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A class of diffeomorphism invariant, physical observables, so-called astrometric observables, is 
introduced. A particularly simple example, the time delay, which expresses the difference between 
two initially synchronized proper time clocks in relative inertial motion, is analyzed in detail. It is 
found to satisfy some interesting inequalities related to the causal structure of classical Lorentzian 
spacetimes. Thus it can serve as a probe of causal structure and, in particular, of violations of 
causality. A quantum model of this observable as well as the calculation of its variance due to 
vacuum fluctuations in quantum linearized gravity are sketched. The question of whether the causal 
inequalities are still satisfied by quantized gravity, which is pertinent to the nature of causality in 
quantum gravity, is raised, but it is shown that perturbative calculations cannot provide a definite 
answer. Some potential applications of astrometric observables in quantum gravity are discussed. 
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I. INTRODUCTION 

The issue of physical observables in both classical 
and quantum theories of gravity has been a topic of 
long standing interest both practically and theoretically. 
Practically, precise models for tracking the positions of 
objects on scales from Solar System to cosmological re- 
quire input from general relativity. Such models go un- 
der the generic name of relativistic astrometry [Tj. Also, 
models of the very early Universe rely on incorporating 
quantum gravitational effects in order to predict poten- 
tially observable signatures in current cosmological ob- 
servations (Sec. 5.7 of [2]). Theoretically, the algebra of 
physical observables, with each observable a mathemati- 
cal model of an experimental outcome, is an integral part 
of a complete classical or quantum theory of gravity [SHS] . 
In each of these cases, one has to confront the problem 
that, while physical observables are expected to be in- 
variant under spacetime diffeomorphisms, in the usual 
formulation of general relativity everything is described 
in terms of tensor fields, which are covariant but not in- 
variant under spacetime diffeomorphisms. 

A resolution of the problem of physical observables 
would then involve two parts. First, one has to ex- 
plicitly describe a sufficiently large class of spacetime 
diffeomorphism-invariant (or diff- invariant) quantities 
expressed in terms of tensor fields. Second, one has to 
identify elements of this class that correspond to out- 
comes of some experiments of interest. The literature on 
this subject is extensive [3HS] (see also references therein), 
but no solution has been entirely successful. To illustrate 
the difficulties, consider the following two examples. A 
simple-to-describe class of diff-invariant quantities con- 
sists of the spacetime integrals of the form 

/ fi9A)Vg, (1) 
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where M is the entire spacetime, / is some smooth space- 
time scalar defined only in terms of the metric g and other 
dynamical fields (f>, and Vg is the metric volume form. Un- 
fortunately, even ignoring the issue of convergence of such 
integrals, this class of diff-invariant quantities is not rich 
enough to describe the outcomes of any experiments that 
we are likely to perform (since such experiments would 
necessarily be localized in a finite region of spacetime). 
The other example is more abstract. Consider (formally) 
the physical phase space of general relativity defined as 
the quotient by spacetime diffeomorphisms of the space 
of solutions of Einstein's equations. Ostensibly, any func- 
tion on this space is a diff-invariant quantity, and hence 
a physical observable. Moreover, all physical observables 
are so captured. However, due to the abstract nature 
of this construction, it is not possible to assign a clear 
physical meaning to any element of this class. There is, 
a priori, no effective way to specify an individual element 
of this class or to carry out practical calculations with it. 

The aim of this work is to take a pragmatic approach 
to the explicit construction of physical observables and 
apply it to more theoretical problems like studying the 
causal structure of quantum gravity. From a theoretical 
point of view, the abstract notion of a physical observ- 
able, sketched in the previous paragraph, as a function on 
the physical phase space is quite satisfactory. The main 
problem remaining is to identify observables of interest 
and given them a physical interpretation in terms of a 
modeled experimental outcome. A natural way of ad- 
dressing this difficulty, inspired by the methods used in 
practical problems like relativistic astrometry, is to start 
with a potential experiment in mind and construct a suf- 
ficiently detailed mathematical model of it. Such a model 
should include sufficiently many dynamical variables rep- 
resenting parts of the experimental apparatus such that 
the desired measurement outcome can be modeled using 
the relative configurations of these variables. The result 
is a mathematical model of a measurement outcome, in 
other words a physical observable. This observable, by 
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virtue of its operational definition, should then be a diff- 
invariant quantity and thus an element of the algebra of 
functions on the physical phase space of the theory. Now 
though, fortunately, since we started out by modeling an 
experiment, its physical interpretation is clear. 

A mathematical model of an experiment interacting 
with dynamical gravity is likely to make reference to so- 
lutions of geodesic or wave equations on unspecified (in- 
deed dynamical) metric backgrounds. Coupled with the 
large variety of experiments that could be imagined and 
modeled, a remaining practical difficulty is that the re- 
sulting physical observable is still specified only implic- 
itly and may not be immediately amenable to practical 
calculations. It appears that this problem can only be 
overcome on a case by case basis. For the particular 
observable considered in this work, the time delay, this 
difficulty is overcome by appealing to perturbation the- 
ory and providing explicit formulas, based on an exact 
implicit definition, in terms of one-dimensional integrals 
over linear metric perturbations about Minkowski space. 

The idea of using operationally or "relationally" de- 
fined observables in gravitational theories is not entirely 
novel. It has been previously considered in [51 [TUl HJ ■ 
Unfortunately, that work has remained at a rather ab- 
stract level and did not make use of sufficiently realistic 
experimental models, thus keeping the physical interpre- 
tation of the constructed observables somewhat moot. 
The previous works that used ideas most similar to ours 
are [12], [13], and [M]. Unfortunately, the original work 
of [15] and its follow-ups [T51 - fT7] , while exhibiting a clear 
physical interpretation, left many mathematical loose 
ends. In particular, the issues of diff-invariance (or gauge 
invariance) and regularization were not treated entirely 
satisfactorily, both of which are explicitly addressed in 
this work, see Sees. VB3 and VI C 2| Another work in 
a similar spirit is jlSj. especially at the technical level, 
though with a different physical motivation. On the other 
hand, the original work |19J gives several different motiva- 
tions for the technical calculations, including an approach 
very similar to that of this paper in terms of the construc- 
tion of diff-invariant, physically meaningful gravitational 
observables. At the technical level, the main departure 
of this work from that of [TH [TH] is in the use of smeared 
observables to regularize divergences appearing due to 
the use of geodesies of idealized, point-like particles, see 

Sec. lvnc] 

In SecjllJ we operationally define the time delay physi- 
cal observable (or rather a family of related observables). 
Section [TIT] gives an exact, though implicit, mathematical 
model for this physical observable in a theory of gravity 
coupled to a minimal amount of matter modeling the 
experimental apparatus. Section |IV| contains an analy- 
sis of why the time delay is an observable interesting for 
studying the causal structure of gravity. In particular, 
two important inequalities are derived directly from the 
Lorentzian character of the metric field. Section jVj us- 
ing technical results on the perturbative solution of the 
geodesic and parallel transport equations presented in 



the [Appendix] gives an explicit formula for the time de- 
lay in linearized gravity. Sections VI and VII sketch how 
the time delay should be defined as a quantum observ- 
able and how explicit calculations in linearized quantum 
gravity can capture some aspects of causal structure of 
quantum gravity. Due to the added complexity of quan- 
tum mechanics, these two sections are naturally less de- 
tailed than the preceding ones. The issues discussed in 
these sections will be addressed in more detail elsewhere. 
Section IVIIII concludes with a discussion of the results 
and an outlook to future work. 



II. OPERATIONAL DESCRIPTION AND 
GAUGE INVARIANCE 

The time delay observable is defined by the following 
experimental protocol, Fig. [T] (which is of course only an 
idealization of a real experiment). Consider a laboratory 
in inertial motion (free fall). The laboratory carries a 
clock that measures the proper time along its trajectory. 
The laboratory also carries an orthogonal frame, which 
is parallel-transported along the lab's worldline. (The 
frame could be Fermi- Walker-transported if the motion 
were not inertial.) At a moment of the experimenter's 
choosing, the lab ejects a probe in a predetermined direc- 
tion, fixed with respect to the lab's orthogonal frame and 
with a predetermined relative velocity. The probe then 
continues to move inertially and carries its own proper 
time clock. The two clocks are synchronized to at the 
ejection event O. After ejection, the probe continuously 
broadcasts its own proper time (time stamped signals), 
in all directions using an electromagnetic signal (which 
hence travels at the speed of light). At a predetermined 
proper time interval s after ejection, event Q, the lab 
records the probe signal and its emission time stamp t(s), 
sent from event P. Call s the reception time, r(s) the 
emission time and the difference 



6t{s) = s — t(s) 



(2) 



the time delay. 

The time delay St{s) as well as the emission and re- 
ception time are presumed to have been measured with 
negligible inaccuracy. Of course, that is a severe ideal- 
ization. For it to be reasonable, the magnitude of St{s) 
must exceed the noise from the intrinsic inaccuracies in 
the probe and lab instruments (clocks, gyroscopes, ejec- 
tion mechanism, transmission and recording uncertain- 
ties, etc.). Given the smallness of both the classical and 
quantum contributions to St{s) (which are suppressed by 
all of the following: magnitude of light speed, smallness 
of spacetime curvature, and smallness of h), is unlikely 
to be reasonable in our own Universe, at least for naive 
ways to realize this experimental setup. 

However, there is no a priori reason for not being able 
to perform such measurements successfully with (signifi- 
cantly) more clever or improved experimental techniques. 
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FIG. 1. Geometry of the experimental protocol defining the 
reception time s, emission time r(s), and the time delay 
St{s) = s — t{s). The synchronization/ejection point is O. 
The signal emission point is P and the signal reception point 
is Q. 



or in a universe with different values of some of the funda- 
mental constants. A successful theory of quantum grav- 
ity should be able to yield quantitative predictions (for 
any universe) for this and related observables. Some set 
of these observables may actually be practically measur- 
able in our own Universe. As such, the time delay, by 
virtue of its simplicity and ease of physical interpretation, 
serves as a useful benchmark for dealing with whatever 
practical difhculties are likely come up in calculations in- 
volving similar, but perhaps more realistic, observables. 



III. CLASSICAL MATHEMATICAL MODEL 

A significantly idealized mathematical model of the ex- 
periment described in the preceding section, in the clas- 
sical theory, consists of the geometrical objects collected 
in the following definition. 

Definition 1. A lab-equipped spacetime {M, g,0,ef) 
consists of an oriented, Lorentzian, time- oriented, glob- 
ally hyperbolic, n-dimensional spacetime {M,g), a point 
O G M and an oriented orthonormal frame ef £ TqM , 
with a an abstract tensor index and i = 0,1, . . . ,n, where 
Cq is timelike and future- directed. 

See Sec. 2.4 of [20] for the distinction between abstract 
and coordinate tensor indices. Physically, the point O 
represents the spacetime event when the probe is ejected 
from the lab. The vector eg is tangent to the lab's world- 
line and ef, i — 1,2, ... ,n is the oriented spatial frame 
carried by the lab. Note that the restrictions on the 
Lorentzian geometry of (M, g) may not all be necessary. 
Also, in this work we only consider the case n = 4. 



Of course, due to the background independence of 
gravitational physics, the measurements carried out in 
spacetimes related by diffeomorphisms (i.e., gauge trans- 
formations) must be identical. It is thus useful to intro- 
duce the following notion of equivalence. 

Definition 2. Two lab-equipped spacetimes (M, g, O, ef) 
and {M' ,g' ,0' ,e'°') are gauge equivalent if there exists 



a diffeomorphism X'- ~^ 
X{0) = O' and x*ef = , 
ferential push-forward. 



M' such that x*9 = ff'; 
where x* denotes the dif- 



An observable is modeled mathematically by a func- 
tion on the space of lab-equipped spacetimes. The time 
delay observable is defined by implementing the proto- 
col outlined in the preceding section. First, we need the 
following further definitions. 

Definition 3. The lab worldline, Q{s) = expQ(su), 
is the geodesic passing through O with tangent vector 

= Cg. Here exp^ : TqM — )■ M is the usual geodesic 
exponential map. 

The probe worldline, P(t) = expo(tu), is the geodesic 
passing through O with tangent vector v"' — v^ef, with 



pi, 3 



a timelike, future- directed, unit vector, chosen 
independent of the spacetime geometry. 

The signal worldline, Z{t,\), is the null geodesic em- 
anating from a point on the probe worldline, Z{t, 1) = 
P{t), and intersecting the lab worldline, Z{t,0) = Q{s), 
with the earliest possible s (alternatively, if s is fixed, 
then t is chosen to be the latest possible). 

The values of t and s connected by Z are functionally 
related. This relationship defines the time delay observ- 
able. 

Definition 4. // t and s are such that P{t) and Q{s) 
are connected by Z(t, X), they are referred to as a pair of 
emission and reception times. The functional relation- 
ship between them is denoted 



t = Ty{s), 



(3) 



where t^(s) is called the recorded emission time and the 
difference 



STy{s) = S - Tv{s), 

is called the time delay. 

By construction, the following theorem holds. 



(4) 



Theorem 1. Given two gauge- equivalent lab-equipped 
spacetimes {M, g,0,ef) and {M' ,g' ,0' ,e'°^), the corre- 
sponding time delays (keeping s and fixed) are equal: 



Sr,{s) = 6t:^{s). 



(5) 



In other words, the time delay (as well as as any func- 
tion thereof, such as the recorded emission time) consti- 
tutes a genuine (diffeomorphism-invariant) physical ob- 
servable on lab-equipped spacetimes. When the context 
is clear, we will omit the explicit dependence of (5t„(s) on 
V or s. 
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IV. CAUSAL INEQUALITIES 

The time delay is interesting in more ways than just 
being an explicit example of a physical observable sensi- 
tive to the ambient gravitational field. The experimen- 
tal protocol defining it can be thought of as designed 
to test the impossibility of superluminal signal propa- 
gation and the geodesic character of the lab worldline. 
In particular, under quite generic assumptions related to 
the Lorentzian character of spacetime geometry, the time 
delay satisfies inequalities that would be violated if the 
above-mentioned tests were to fail. A careful examina- 
tion of a mathematical model of such tests for classical 
spacetimes can serve as a benchmark to understand pos- 
sible outcomes of such tests in any proposed theory of 
quantum spacetime. 

A. Maximality of light speed 

In curved spacetimes, it is impossible to objectively 
compare the speed of light at different spacetime events. 
Due to general covariance, any experiment to measure 
the local speed of light, calibrated at spacetime event x 
to return the value 1, will return the same value 1 at 
any other spacetime event y, provided it was parallel- 
transported there. Therefore any such experiment must 
perform measurements in finite regions of spacetime. 

In the case of the time delay experiment, since light 
is used to send the signals, the local speed of signal 
propagation cannot by definition exceed that of light. 
However, to take into account possible nontrivial global 
geometry, we adopt the following definition of (appar- 
ent) superluminal signal propagation. Consider a pair 
of emission-reception times (r, s) . If there exists another 
pair (t', s') such that t' < t (signal emitted later than 
t') and s' > s (signal arrived earlier than s'), then the 
later signal must have travelled superluminally, cf. Fig. [2] 
In classical Lorentzian geometries without closed causal 
curves we can prove that this never happens. 

Theorem 2. In a lab-equipped spacetime {M, g^O,ef) 
we have the following implication between inequalities 
satisfied by pairs (r(s),s) and (t(s'),s') of emission- 
reception times: 

t(s') < t(s) =^ s' < s. (6) 

In particular, when t(s) is smooth, we have ^t{s) > 0. 

Essentially, this theorem says that a signal that is emit- 
ted later, with respect to the probe, also arrives later, 
with respect to the lab. 

Proof. Let the two signals be emitted from chronologi- 
cally successive points P' and P, and received at points 
Q' and Q, respectively. Since P' and P are part of the 
same worldline, P clearly belongs to the set of all points 
that can be reached from P' by future-directed timelike 



Q 
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FIG. 2. Illustration of the conclusion of Theorem |2l Succes- 
sive emission times t < t imply successive reception times 
s < s' . The dashed line represents a case, ruled out by the 
theorem, where the signal might appear superluminal. 

curves, P G /+(P'). The points Q and Q' are also con- 
nected by a timelike curve, though we do not assume in 
which precedence order, therefore Q must belong to ei- 
ther I^{Q') or I^{Q'). At the same time, by the defini- 
tion of Q, it can be reached by a piecewise smooth, non- 
spacelike, future-directed curve P'PQ. Since P'PQ is 
obviously not a null geodesic. Prop. 4.5.10 of [21] implies 
that Q and P' can be joined by a (future-directed) time- 
like curve, Q e /+(P') C int J+{P'). By definition, Q' is 
reached from P' by a future-directed null geodesic, such 
that there is no later point P" on the probe worldline 
with the same property. This implies that Q' S dJ'^{P'). 
Otherwise, Q' G /+(P'): hence P' S I~{Q'), hence any 
point P" S dJ~{Q') that is also on the probe worldline 
violates the preceding hypothesis. But all the timelike 
curves from 9J+(P') to int J+(P') 3 Q can only be fu- 
ture directed, of which the one reaching Q from Q' is 
a special case, hence Q G I^{Q'). This shows that Q' 
chronologically precedes Q oi s' < s. □ 



B. Geodesic extremality 

The twin "paradox" is a well-known phenomenon in 
special relativity: the proper time between two timelike 
separated events is maximized by a straight line (inertial 
motion). Its generalization to curved spacetime is gen- 
erally true only locally: a timelike geodesic maximizes 
proper time among causal curves close to it (provided it 
has no conjugate points). Under some conditions on the 
spacetime or under some extra restrictions on the class 
of allowed causal curves, geodesic extremality can also 
hold globally. This includes the special geometry of the 
time delay experiment. As we shall see below, the time 
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FIG. 3. Illustration of the proof of Theorem [3] The auxiliary 
dashed curve interpolates between OQ and OP, as A varies 
from to 1. As it does so, its proper time length T(A) is 
shown to decrease monotonically, thus implying t(s) < s. 



delay observable is also sensitive to some violations of 
the geodesic extremality. Such violations mimic a break- 
down of the equivalence principle (objects no longer fall 
on geodesies in the absence of external forces). 

Theorem 3. In a lab-equipped spacetime {M, g,0,ef) 
(where the lab and probe worldlines are smoothly de- 
formable into each other in a sense to be precised in the 
proof) a pair (r(s),s) of emission-reception times satis- 
fies the inequality 



t{s) < s or equivalently (5t(s) < 0. 



(7) 



Proof. The basic strategy of the proof is to construct a 
one-parameter family of piecewise geodesic curves that 
interpolate between the lab worldline OQ and the probe- 
signal worldline OPQ, while their proper time lengths 
decrease monotonically, cf. Fig. [3] The existence of the 
specific interpolation constructed below is the extra tech- 
nical hypothesis alluded to in the statement of the theo- 
rem. 

Suppose that the PQ geodesic is affinely parametrized 
as z\\), where Z(0) = Q and Z{1) = P. Denote also 
Pa = Z{X). Then the family OP\Q clearly interpolates 
between OQ and OPQ, where OPx is a timelike geodesic 
connecting these points and P\Q is the segment Z{[Q, A]). 
Since the P\Q segment is null, only the OP\ segment 
contributes to the proper time T{\) along OP\Q. Since 
T(0) = s and r(l) = t(s), the proof is concluded as soon 
as we show that g^T(A) < 0, which we do below. 

We adapt the calculation of the first variation of the 
proper time length of a piecewise geodesic curve from 
Prop. 4.5.4 of |2l]. Let Y{t,\) denote the geodesic 
family OP\, parametrized such that y(0. A) = O and 
A) = Z{\). Y{t, A) is assumed to be smooth in both 



arguments by the smooth deformability hypothesis of the 
theorem. Denote = §-^Y{t,\), Y'" = ^Y(t,X) and 

/ = [-5„fc(r(t,A))y°rY^^- Also, for the purposes of 
the calculation below, pick a coordinate chart a;" and 
replace the Latin abstract tensor indices by Greek coor- 
dinate indices. 



/' = -r 



yf^ 1 



Y'iy°'Y^ 



d 



Y0 



(8) 
(9) 
(10) 




1 d 



YP 

-fdt l-^-^^T 



1 Y°'Y^ 



dt \ ^""-^'^ / 



(11) 



Note that the bracketed term vanished because it is 
precisely the geodesic condition (Eq. 87.3a of [22j) and 
Y{t, A) is a geodesic for fixed A. At t = 0, we have 
Y'" = 0, while at t = 1, we have Y'" = which is a 
past-directed null vector. 



T'(A) 



d_ 
dX 



r dti-gabY'^Y"]'^^ ^ f dtf (12) 

(13) 



^ d Y^ 



dt 



Yb 

-gabZ^'—KO. 



(14) 



The latter inequality follows because Y°' j f is a future- 
directed timelike unit vector and Z'" is a past-directed 
null vector, hence their inner product is positive. Armed 
with this inequality, it immediately follows that 

s - t{s) = T(0) - T(l) = - / dA T'(A) > 0, (15) 

Jo 

which completes the proof. □ 



V. EXPLICIT CALCULATION IN CLASSICAL 
LINEARIZED GRAVITY 

The time delay observable, while well-defined from 
its description in the preceding sections, has so far 
been defined only implicitly. Unfortunately, it would be 
very difficult to obtain an explicit expression for it, ex- 
cept in highly symmetric spacetimes, where the required 
geodesies can be computed explicitly. In particular, in 
Minkowski space, as is done below, it can be computed by 
elementary means. Fortunately, for small perturbations 
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of Minkowski space, an explicit expression for the time 
delay can be found at linear order. Such an expression 
would be especially needed for the calculation of quan- 
tum averages and fluctuations, as sketched in Sec. |VII| 

The calculations are carried out in the tetrad formal- 
ism. While linearized gravity calculations are usually car- 
ried out in the more familiar metric variables, there are a 
few reasons to consider tetrads. Using tetrads opens the 
door to a kind of improved perturbation theory, where 
the metric keeps its Lorentzian signature at every step of 
the approximation. This line of investigation, as briefly 
brought up in Sec. |VIII[ will be pursued elsewhere. An- 
other advantage of tetrads is that they are needed in the 
standard way of formulating fermions on curved space- 
time. 

First, we explicitly compute the time delay in Minkow- 
ski space and check the causal inequalities. Then, using 
the results of the perturbative solution of the geodesic 
and parallel transport equations of the [Appendix} we 
compute the explicit expression for the time delay at lin- 
ear order in the deviation from Minkowski space. 



Note that 9 > since the probe is moving away from 
the lab. The null vector connecting the emission and 
absorption points is 

w' = su' - tv' = s{u' - e-\'). (23) 

Another useful identity is 

v' = e\u' -w^js). (24) 

B. Approximately Minkowski space 

1. Tetrad formalism 

Consider another lab-equipped spacetime (M, g, 0, ef), 
where we have kept the same underlying manifold M 
and synchronization point O = as in Minkowski space. 
On the other hand, we express the new metric as gab = 
Vij^a^i^ where and e° is a new dual pair of orthonormal 
tetrads. 



A. Minkowski space 

Consider Minkowski space (M = K^, 77, 0, if ), with 
r] = diag(— 1, 1, 1, 1), as a lab-equipped spacetime. With- 
out loss of generality, we can take an arbitrary inertial 
coordinates on (ill, 77) and use their origin as the 
synchronization point and the vectors xf = {d/dx^)°' as 
the reference tetrad. The dual tetrad is x^ — (da;')a and 



satisfies the identities xfa 



Si and x'^xl ^ 5l. 



The 



Minkowski metric is rjab = Vij^a^i- 

The lab and probe worldlines are parametrized, respec- 
tively, as x'(s) = su* and x^{t) = tv^, i = 0,1, 2, 3, where 

= (1,0,0,0). Suppose that the relative speed of the 
two timelike vectors it" = u^x"; and v"- — v'^xf is given by 
the positive hyperbolic rapidity 9, Wroi — tanh^, then we 
have u - V = rjabu^'v^ = "nijU^v^ = — cosh0. The values 
of s and t which may be connected by light signals are 
constrained by 



riij{su^ — tv^){su^ ~ tv^) — 
-s^ - 2st{u ■v)-t^ = 
- 2st cosh 6* + = 
(se^ -t)(se"'' -i) = 0. 
The retarded solution is then 

Tci(s) = t = se~^. 



(16) 
(17) 
(18) 
(19) 

(20) 



The subscript stands for "classical," as it will serve in 
Sec. |VII| as the classical background expectation for quan- 
tum fluctuations. This expression clearly satisfies the 
causal inequalities obtained in the previous section: 



ds 



Tcl(s) 
Tcl(s) 



se 



" < s, 
> 0. 



(21) 
(22) 



(25) 



Using the Minkowski tetrad on M as a reference, any 
other one can be obtained by a local general linear trans- 
formation 



a _ rpi -o i _ rpi -i' 



(26) 



where T and T are spacetime-dependent invertible ma- 
trices, such that T = T~^. Similarly, any lab frame e° 
can be obtained by another general linear transformation 
at O, 



e^ = {To)lxf,. 



(27) 



The possible discrepancy between the lab frame and the 
spacetime tetrad at O is 



(28) 



where L is clearly a Lorentz transformation, L\ 7]iiji — 
Vij- 

If this new lab-equipped spacetime is approximately 
Minkowski, then both T^, and must be close to the 
identity matrix. This is conveniently expressed by first 
parametrizing them as T = exp(/i) and L = exp(/io), 
and then requiring that h and ho are close to 0. The 
smallness requirement aside, h and ho could be, respec- 
tively, an arbitrary matrix and an arbitrary skew-adjoint 
matrix, ?7ifc(/io)j = ~Vkj{ho)i ■ Then the metric is 



r]ij e j Bj 



n,,Tl,Tl,x^xl (29) 

= lab + [v^■3K' + mr^'WlH + o{h^) (30) 

= rjab + hab- (31) 
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The last two equations describe the relationship between 
the deviations /i* and hab from Minkowski space, in the 
tetrad and metric formalisms respectively, 

hab — '2'h(ij)X^a^l + 0{h?), where hi 



(32) 



A worldhne ^{t) is described by its coordinates 7*(t) = 
x^{"f{t)). Its tangent vector is denoted 7(f)°. Knowledge 
of the tangent vector allows one to recover the curve as 
follows 

\tr{t){dx')a= r^'\x' =^\h)-i\h). (33) 

L "'7(*l) 

For convenience, all curves are afhnely parametrized from 
to 1. Thus, the length of a timelike geodesic is equal 
to the length of its initial tangent vector. 

A geodesic ^{t) is completely specified by its point of 
origin 7(0) and its initial tangent vector 7''(0), while a 7- 
parallel-transported vector v°'{t) is specified by its initial 
value ^"(O) at 7(0). Again, for convenience in further cal- 
culations, all such initial data are specified with reference 
to some given curve /3, with /3(0) = O. Namely, the point 
of origin is 7(0) = /?(!), the initial tangent vector 7°(0) is 
the /3-parallel-transported image of a vector 7q = jp^f^ 
and the initial value ^"(O) is the /3-parallel-transported 
image of a vector Vq = Vpef (cf. Fig. [7]). The geodesic 
and parallel transport equations are written down and 
solved to order 0{h) in the Appendix 



triangle closes (the end point of U is in fact O), and fi- 
nally that the tangent to C/ at O is —su^ef (which is also 
the parallel-transported image along the VWU triangle, 
in other words a holonomy image, of —su°'): 

S = s, (34) 

i = e^t, (35) 

= „*e° =e,"exp(;^o)>^ (36) 

u° = e^[exp(p,7) exp{pw) exp(p^/)]j(/^o)i^i^ (37) 

u;" = e^exp(g)>^ (38) 

where we have used the notatio n exp( p^) for the parallel 
transport operator along 7, Eq. (A. 24), while f is a scalar 



and exp(q') a Lorentz transformation (rjikq 
both yet to be determined. Note that does not 
parametrize w"' uniquely, as exp(q) could always be pre- 
multiplied by another Lorentz transformation fixing w'', 
but it does contain three non-arbitrary parameters. The 
only condition left to be satisfied is the closure of the 
VWU triangle (equating the end point of U with O), 
which provides four equations. These four equations can 
be used to solve for the remaining undetermined param- 
eters, one in f and three in gj. Since we are working 
at linear order, we only need the leading terms in the 
expansion of these unknowns 



9, = 9, 



0{h'), 



0{h^ 



(39) 



2. Geodestc tnangle construction 

All curves considered in this section are perturbations 
of piecewise linear paths, which are piecewise geodesic in 
Minkowski space. In particular, at zeroth order in h, the 
sides of the geodesic triangle formed by the worldlines 
of the lab, the probe, and the signal form an ordered 
sequence of spacetime segments {V,W,U), as illustrated 
in Fig. [4j Namely, V stretches from O to P,W stretches 
from P to Q, and U stretches from Q back to O. Using 
the convention of the last paragraph of the preceding 
section, each of the {V, W, U) segments can be specified 
as starting from the end point of the preceding one (note 
that the order corresponds to counterclockwise starting 
from O in Fig. [4| with the respective tangent vectors 
{tv"" ,w°' ,—su°'). Because Minkowski space is flat, it is 
clear that the segments VWU form a closed triangle by 
virtue of their tangent vectors adding up to zero. 

In approximately Minkowski space, we wish to describe 
a perturbed version of the above construction. Namely, a 
sequence of geodesic segments {V, W, U), connected from 
end to end, with the respective images {tv°',w°',—su°') 
of their initial tangent vectors parallel-transported to O. 
We take and v"" to be unit vectors, hence s and i 
are the proper time lengths of the corresponding seg- 
ments. To be consistent with the experimental protocol 
described in Sees. |TT] and |III[ we must take s = s and 
V = v^ef, require that li'^ is null, require that the geodesic 



Using the perturbative solution of the geodesic and 
parallel transport equations obtained in the Appendix 
(Eqs. (A.30) and (A.34)), at linear order, the triangle 



closure condition can be written out explicitly as 
+ («/ + q^w^ + J\vy) 



(40) 



(41) 



where, using the notation of Eqs. (A.57) and (A.58), we 
have defined 



VikH^ — Hij — iH(^v,w,u))ij7 

VikJ^ ^ Ji — iJ(V,W,U))i- 



(42) 
(43) 



The expression in parentheses vanishes due to the closure 
of the zeroth-order geodesic triangle. Also, contracting 
the closure condition with w' makes the term with q van- 
ish (due to its antisymmetry). The solution for r is then 



ui'Ji -I- w^HijSU^ 

Tclis)v ■ W 



(44) 



The detailed structure of the defining expression for 
r in Eq. ( 44 ) can be deduced from the structure of the 



expressions for the H and J terms, given explicitly in 
Eqs. (A. 43 1 and (A.46). It can be described as follows. 



FIG. 4. Schematic structure of the ff-term in r, Eq. (46 1. 
Notation foUows Eqs. ((46|-([49|. 



FIG. 5. Schematic structure of the J-term in r, Eq. (47 1. 
Notation follows Eqs. ((46|-(|49|. 



Both H and J consist of a sum of terms associated to 
the segments of the VWU triangle. A term associated to 
segment X consists of a tensor, built up from the vectors 
u*, and w% contracted with a (possibly iterated) line 
integral over X, where the integrand consists of the per- 
turbation hij, possibly with several derivatives applied 
to it. Schematically, this structure can be expressed as 



n(m) 



X 



(45) 



X 



where all tensor indices have are suppressed and iter- 
ated integrals are represented using the notation from 
Eqs. ( A.37)-( A.41 ). There is at most one derivative 
{k < 1) and integration over a spacetime segment is iter- 
ated at most twice (m < 1). In a bit more detail, though 
leaving the tensor contractions aside, the structure of the 
H and J terms can be expressed as follows 



H 



x=v,u,w 




Y<X 



(46) 

(47) 

(48) 
(49) 



where the order between the segments is counterclockwise 
starting from O, as in Figs. |4] and [5| The geometry of 
the various terms is illustrated in Figs. |4] and [5] This 
information is used in Sec. IVII Dl 



3. Time delay and gauge mvariance 

As proven in Theorem [l] the time delay is a gauge- 
invariant observable. From the formula 

t(s) = i = eVei(s) = Tci(s) +rrei(s) +0(/i2), (50) 

that relates the linearized gravity correction rT(s), 



Eq. (44), to the Minkowski space result rci(s), Eq. (20 1, 
it is obvious that r should be invariant under linearized 
gauge transformations. This can be checked explicitly 
using the gauge transformation formulas, Eqs. (A.52) 



and (A. 56 1, for the terms making up Hij and Ji. As 



a consequence, which is given at the bottom of the |Ap- 
|pendix[ the closure of the VWU triangle in Minkowski 
space implies the individual gauge invariances of both 
Hij and Ji, and hence of r. 

The last remark deserves some emphasis. There have 
been many attempts to try to achieve some sort of ex- 
plicit and complete classification of gauge-invariant ob- 
servables of general relativity |3H11 [!]• So far, no such 
complete classification is known. Even in the case of a 
partial classification, such lists of gauge-invariant observ- 
ables are often obtained without direct physical interpre- 
tation. The strategy of this paper has been different. 
The idea was to first establish an operational definition 
of an observable, in terms of the thought experiment de- 
scribed in Sec. |Tlj second to establish a mathematical 
model thereof, which would naturally be gauge-invariant 
though perhaps only defined implicitly, and third to use 
an approximation method (linear-order perturbation the- 
ory, in this case) to obtain an explicit expression for the 
observable. The result of this strategy is an explicit (lin- 
early) gauge-invariant expression for an observable r and 
a physical interpretation of it as an approximation to the 
outcome of a clearly described thought experiment. It is 
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of course highly hkely that an exhaustive classification of 
gauge-invariant observables, for the simpler problem of 
linearized gravity, would have identified explicit expres- 
sions like Hij and Ji, but it is at the same time highly 
doubtful that they would be accompanied by the clear 
physical interpretation we have managed to associated 
to their particular combination in ( 44 ) . 

It is also worth noting that the works of Ford et 
al. [m [T5HT7] and Roura and Arteaga [13] worked in 
a particular gauge and with more restricted experiment 
geometries. Thus they did not obtain the same gen- 
eral gauge-invariant expressions that we have derived 
here. However, similar expressions, expanded even to 
quadratic order, were obtained in the work of Tsamis 
and Woodard [T8l. 



VI. SKETCH OF QUANTUM MATHEMATICAL 
MODEL 

Ideally, to be able to theoretically describe quantum 
effects, the thought experiment protocol described in 
Sec. HI] should be translated into a mathematical model 
within a quantum theory that encompasses both the 
gravitational field and the experimental apparatus de- 
scribed in the protocol. A naive attempt to do this is 
obstructed by several difficulties: (a) the lack of a uni- 
formly accepted (or at the very least sufficiently general) 
quantum theory of gravity, (b) the identification of a time 
observable in quantum mechanics, and (c) the difficulties 
in modeling measurements in quantum mechanics. For- 
tunately, we can propose pragmatic solutions to each of 
these problems, as discussed below. 



A. Quantum linearized gravity 

While it is true that there is no uniformly accepted 
theory of quantum gravity, there are some common stan- 
dards that are expected to be met by the final version of 
any proposal. One such routine benchmark is the ability 
to reproduce classical general relativity in the appropri- 
ate limit. It is worth noting that under very general 
circumstances (in the absence of strong curvatures), the 
dynamics of the gravitational field in general relativity 
can be very closely approximated by the dynamics of lin- 
earized gravity, also known as the theory of (linear) grav- 
itational waves. Our experience to date overwhelmingly 
demonstrates that the quantum theory of any field whose 
dynamics may be approximated by a linear theory, be it 
a "fundamental" field as in elementary particle physics or 
an "effective" field as in condensed matter theory, is well- 
approximated by the Fock quantization of the approxi- 
mate linear theory. By inductive reasoning, we presume 
that any proposed theory of quantum gravity should also 
be benchmarked by its ability to reproduce quantum lin- 
earized gravity. Therefore, pragmatically, we restrict our- 
selves to the Fock quantization of the linearized gravity 



field on Minkowski space as the approximate quantum 
theory of gravity for the purposes of the mathematical 
model of the time delay observable. 



B. Time in quantum meclianics 

It is often repeated physics lore that there is no ob- 
servable in quantum mechanics corresponding to time, 
which naturally leads one to wonder whether it is even 
possible to model time measurements in quantum me- 
chanics. This argument is originally due to Pauli (p. 63, 
footnote 2 of [23 )• Fortunately, when precisely stated, 
it is much less restrictive than one is first lead to be- 
lieve pH - E7] . The crux of this argument is a contradic- 
tion that stems from the following hypotheses. Suppose 
we have a quantum mechanical system with Hamiltonian 
H, whose spectrum is bounded from below, and an oper- 
ator observable T, whose commutation relation with H is 
precisely of the form [T, H] = ih, as would be appropri- 
ate for a "time observable" T (together with appropriate 
continuity and functional analytical conditions). Then, 
an appeal to the Stone- von Neumann uniqueness theorem 
(Theorem VIII. 14 in [3S]) establishes a contradiction, as, 
according to the theorem, both T and H must have con- 
tinuous unbounded spectra. Thus, there cannot exist 
such an observable T corresponding to time. However, 
there are at least two physically reasonable ways to cir- 
cumvent this conclusion. One is to drop the hypothesis 
that H is bounded from below. While this requirement is 
important for the global, long-term stability of physical 
systems, its not necessary in some approximate descrip- 
tions meant to describe the dynamics of some system 
for bounded time intervals. Two common examples are 
a particle in a linear potential and a harmonic oscilla- 
tor with an inverted potential. The other is to relax the 
commutation relation condition to [T, H] « ih, where the 
correction terms that restore equality may be higher or- 
der in h or may be small in another way when restricted 
to a physically relevant subspace of possible states. An 
example is a particle on a circle, whose dynamics dictate 
uniform motion, so that its position can serve as an ap- 
proximate "cyclic time" observable, like the position of 
the hand of an analog clock. Many more examples are 
discussed in and the references therein. 



C. Modeling quantum measurements 

1. Classical vs quantum measurements 

The remaining obstacle is overcome by constructing a 
fairly explicit, though still rough, model of a measure- 
ment, where the system of interest (gravitational field, 
lab, probe, signal), the measurement devices (proper time 
clocks) and recording devices are all taken into account. 
The details of this setup are described below, following 
some of the ideas of j^H 1301 on the use of physical clocks 
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in quantum systems. The conclusion can be formulated 
as follows. After the reception and emission times of 
a signal have been measured by the lab and individually 
stored, the states and the dynamics of the storage devices 
stabilize and decouple from the rest of the system, as well 
as from each other, in the asymptotic future. Then, in 
the asymptotic future, the corresponding "readout" ob- 
servables S (recorded reception time) and T (recorded 
emission time) commute and thus define a joint (classi- 
cal) probability distribution p{a,t) = {6{S — (t)6{T — t)) , 
where the expectation value is taken with respect to the 
(Heisenberg) state of the total system, which we will refer 
to as the quantum gravitational vacuum. Mathematically, 
this probability distribution p{a, t) may be referred to as 
either the joint spectral density of the quantum gravita- 
tional vacuum with respect to the operators S and T. In 
more physical terms, p{a, t) is the absolute value squared 
of the wave function of the quantum gravitational vac- 
uum projected onto the variables a and t. 

Recall that the main output the classical mathematical 
model of the measurement of the time delay observable 
is the functional relation t = t{s), which can be seen as a 
special case of a joint probability distribution pc\{a,t) — 
5{a — s)S{t — t((t)), where s is the predetermined time 
when the laboratory makes the measurements. This clas- 
sical probability distribution is so "sharp" because we 
take as a classical state a definite configuration of the 
gravitational field. More generally, in the framework of 
classical statistical mechanics, we can take any proba- 
bility measure dp{g) on the space of gauge equivalence 
classes of the configurations of lab-equipped spacetimes. 
The main output of the classical mathematical model of 
the time delay measurement in this state is then the prob- 
ability distribution Pel (f, i) = j dp{g)6{cf—s)5{t—Tg{a)), 
where the dependence of the emission time Tg{s) on the 
equivalence class g of lab-equipped spacetime configura- 
tions is indicated through a subscript. Thus, considering 
quantum mechanics as an extension (or rather deforma- 
tion) of classical statistical mechanics, it is not surpris- 
ing that the main output of a quantum mathematical 
model of a measurement of the time delay observable is 
the probability distribution p{(j,t). Of course, being the 
result of a quantum measurement, the distribution p{a, t) 
depends on more details of the measurement (such as the 
order in which the measurements were carried out) than 
the classical distribution Pci{cr, t). 

2. Dynamical apparatus model 

The full system included in the model consists of the 
following dynamical subsystems: the gravitational field 
g, the lab and probe worldline coordinates yi and jjp, the 
lab and probe proper time clocks f; and fp, the time reg- 
isters S and T in the lab, the coordinates of the signal 
particles z, and the time stamp f carried by each signal 
particle. The spacetime is presumed to have a fixed foli- 
ation by level sets of a time function t. The gravitational 



field is taken to be completely gauge fixed, for instance 
using the transverse, traceless, and i-compatible radia- 
tion conditions (Sec. 4.4b of [20] )■ AH worldlines can 
then be parametrized by t as well. The dynamics of the 
full system, describing its evolution with respect to time 
t, is specified by a Hamiltonian 

H = Hsuh + -f^gcom + ^mcasj (^1) 

which is composed of Hguh describing the independent 
dynamics of the subsystems, of Hgeom describing the nec- 
essary interactions or external interventions to effect the 
geometry of the experimental setup, and of -ffmcas de- 
scribing the coupling between the recording devices and 
the rest of the system during the measurement. 

Since we are mostly concerned here with a quantum 
model of measurement of the clock readings, we will 
concentrate only on /fmeas and specify Hsuh and -ffgoom 
mostly verbally. 

The dynamics of the gravitational field follow the ap- 
propriate gauge fixed Hamiltonian, a term in Hguh de- 
rived from the Einstein-Hilbert action. The precise de- 
tails of the implementation of this idea are irrelevant for 
this discussion, as long as the corresponding dynamics 
about the quantum gravitational vacuum can be approx- 
imated by the dynamics of linearized gravity about the 
Fock vacuum. This assumption is the basis of the calcu- 
lation sketched in Sec. IVIII 

The worldlines of various particles are described by 
their spatial coordinates as functions of the global time 
t, yl{t), yp{t), and z^{t). The dynamics of these variables 

follow from the appropriate terms in Hsuh- The lab and 
probe worldlines are timelike geodesies, with an appro- 
priate term in /fgcom providing a kick to the probe at 
event O to give it a fixed relative velocity with respect 
to the lab. (The fact that O lies on the lab worldline 
can be used as one of the gauge-fixing conditions.) To 
imitate the action of a continuously emitted signal field 
(like the electromagnetic field), the multiplicity of signal 
particles are indexed by a time t' and a unit 3-vector n. 
The dynamics as specified by terms in -ffgub and i?geom 
should be as follows. The worldline z^^, ^-^{t) follows the 
probe worldline until the time t — t' , after which point 
the worldline of z^^, (t) becomes null with direction de- 
termined by n. (This is a kind of eikonal approximation, 
which replaces a massless field by a large collection of 
massless particles.) The initial state of each of these par- 
ticles is presumed to be of a localized wave packet form, 
with negligible wave packet spread on time scales com- 
parable to the geometry of the experiment. 

There are two potential problems in constructing a de- 
tailed quantum model implementing the above require- 
ments. While it is not difficult to write down a classical 
version of such -ffsub + ^^gcom , the generalization to quan- 
tum mechanics is not unique, due to the usual operator 
ordering ambiguities. The standard solution of this prob- 
lem is to parametrize these ambiguities and realize that 
different choices of these parameters correspond to physi- 
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cally different models. Thus, the fixation of these param- 
eters must be part of the full specification of the detailed 
model. Fortunately, these ordering ambiguities are gener- 
ically expected to be suppressed by powers of h. More- 
over, we assume that their parametrization may be tuned 
to maximize the validity of the approximations used in 
Sec. |VIl[ The second problem is that coupling point parti- 
cles to fields generically leads to singular dynamics. (The 
singularities inherent in the naive interaction of a classi- 
cal point electron with its own electromagnetic field is a 
classical example of this difficulty.) However, this issue 
can be dealt with straightforwardly by spatial smearing 
of the particle-metric field interaction terms. The spa- 
tial extent of the smearing becomes another parameter 
whose value is to be chosen as to minimize the impact of 
the smearing on the rest of the discussion. Alternatively, 
the interaction term could be modified in a more sophis- 
ticated way, without introducing non-local smearing, for 
instance along the lines suggested by the recent work on 
classical point particles coupled to their self-force [3T] or 
by appealing to intrinsic quantum uncertainty of the cen- 
ter of mass coordinates as in [32] ■ 

The state spaces for the clock and time register sub- 
systems can be presumed to be completely internal (i.e., 
divorced from spacetime coordinates) and thus can be 
subject to even further simplifications. The time register 
subsystems S, T and „) should be very stable, thus 
their contribution to -ffsub should be approximately zero. 
On the other hand, the clock variables fi{t) and Tp{t) 
should evolve approximately monotonically, with rates 
set by their local proper time. This can be accomplished 
by a contribution to i/gub of the form f/P,-, -\-tpPr^ , where 
Pti and Ptj, are respectively canonically conjugate to f; 
and fp, while f/ and fp stand for the appropriate expres- 
sions in terms of y\, {jp and g. Note that this choice of 
Hamiltonian circumvents Pauli's impossibility argument 
by virtue of being unbounded from below. 

Finally, /fmcas is chosen to implement the idea of weak 
measurement [33 . The idea of weak measurement can be 
described as follows. Suppose there is a quantum vari- 
able q whose value we wish to measure and record in 
another variable Q, belonging to a recording device sub- 
system. Suppose that P is canonically conjugate to Q 
and that Q suffers negligible evolution on its own. Then 
the value of Q can be measured at any convenient time 
after the weak measurement took place, thus allowing us 
to infer (subject to quantum uncertainties) the value of q 
at the time of measurement. The measurement itself can 
be modeled using the interaction Hamiltonian ftng{t)qP, 
where ftngit) is a trigger factor, which is non-zero only 
during the time interval when the measurement is sup- 
posed to take place. The operators q{t) and P{t) are 
presumed to commute at equal times, as they belong to 
independent subsystems, so their ordering of the inter- 
action Hamiltonian is unambiguous. If this interval is of 
length At and during it /trig(i) ~ fo is approximately 
constant, we can see that this interaction Hamiltonian 



effects the evolution 



Quite 



- „»/dt/tHg(t)9P 



^before ^ 



Qbcforc + J <itfuigit)q 
Qboforc + Atfoq{to), 



(52) 
(53) 
(54) 



where the last approximation holds provided to was part 
of the measurement time interval and q{t) and P{t) 
evolved negligibly during it. In general, the trigger /trig 
need not be a scalar, and may itself be a operator that 
commutes with both q{t) and P{t). Also while each 
pair of factors commutes at equal times, in general, they 
will not commute at unequal times (even with them- 
selves). Thus, the evolution effected by the interaction 
will involve the time-ordered exponential of the interac- 
tion Hamiltonian and will look more complicated than 



Eq. ( 52 ) . However if the interaction Hamiltonian can be 



considered as a small perturbation, then at linear order 
Qaftcr will look the same as Eq. ( 54 ) . 

With the above discussion in mind, we set the mea- 
surement interaction Hamiltonian to 



-ffmoas ^ S[ti- s]ti 

TiPs+ I dt'dnS^^ 



(55) 



Note that the factors in each product commute at equal 
times (recall that f; is not canonically conjugate to f/), 
so their ordering in iJmcas is unambiguous. The extra 
factor of f;(<) is there to ensure that i?mcas is defined in- 
dependent of the choice of the background time t. Under 
the hypotheses explained in the previous paragraph, the 
asymptotic future values of S and T operators can be 
approximated as 



5+ = lim S{t) 

t—^oo 



(56) 



« ^before + J dtS[Ti{t) ~ s]fi{t)Ti{t) (57) 



r+ = lim T{t) 



Tbeforo + J dt6[Ti{t) - s]f;(t), 

X / dt'dnS^[z^t',r.){t) - yi{t)]rit',r.){t). 



(59) 



Provided S'bcforo and Tbcforc have zero expectation value, 
the measurements of and the above asymptotic 
limits provide unbiased estimates of the remaining terms, 
which can be interpreted respectively as the reception 
and emission times defined in the time delay experimen- 
tal protocol. 

We conclude this analysis by noting that, provided that 
any potential uncertainties can be neglected or modeled 
and subtracted, it is reasonable to assume that the spec- 
tral density of the quantum gravitational vacuum with 
respect to S can be well-approximated by 



{S{S+ - a)) ^ S{<T - s). 



(60) 
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It follows that it is then also reasonable to assume 
that the joint spectral density of the quantum gravita- 
tional vacuum with respect to and T"*" will be well- 
approximated by 

p{<j, t) - {S{S+ - a)5{f+ - t)) « 5{a - .s)p,(t). (61) 

The probability distribution Ps{t) then has the interpre- 
tation of the spectral density of the quantum gravita- 
tional vacuum with respect to the quantum emission time 
operator observable f(s), where the quantum emission 
time observable can then be identified as f (s) = . The 
quantum time delay observable is then simply 5t{s) = 
s — f(s). 

At this point it is worth considering a bit more pre- 
cisely how the probability distributions of Eqs. (60) 



Pintle — s). The joint probability distribution (61 1 would 



and (61) relate to those that would be obtained by a 



physically realized experiment following the same oper- 
ational protocol. Above, we have explicitly stated that 
these expressions are expected to be applicable provided 
all sources of quantum (or classical) fluctuations other 
than the quantum gravitational vacuum are neglected. 
How can this neglect be reasonable? These neglected 
sources are numerous, as for instance discussed in the 



later Sec. VII C Moreover, by now, the study of uncer- 
tainties in measurements of space and time intervals in- 
duced by quantum fluctuations of the internal states of 
a measurement apparatus is classical subject, going back 
to a seminal paper of Salecker and Wigner [34J. These 
effects can be quit e large compared to the Planck-scale 
effects (Sec. VII D ) that we are concerned with here. 

The main difference between the effects we neglect and 
the effect that we actually study is that the former de- 
pend primarily on the internal physics of the appara- 
tus, while the latter depends crucially on the dynamical 
quantum gravitational field. That is, the effect that we 
study is genuinely due to quantum gravity, while those 
we neglect are not. If we are concerned with a question 
of principle, which is to account for all possible sources 
contributing to the variance of the time delay or time 
emission observables described earlier, it is not sufficient 
to include only the internal apparatus sources or only the 
quantum gravitational effects (independent of their rela- 
tive size) , it is in fact necessary to consider both of them. 
The internal apparatus fluctuation sources have already 
been studied extensively in the literature spawned by the 
original work '34]. On the other hand, genuine quantum 
gravitational effects have received much less attention (a 
review of the relevant literature was given in the |Intro- 



duction) and are hence the main focus of the current 
work. Since these contributions to the observational vari- 
ance are separate, they can be analyzed separately and, 
at the leading perturbative order, contribute essentially 
additively. 

When restored, the effects of quantum fluctuations of 
the internal dynamics of the clocks and recording de- 
vices used in described measurement models replace the 



sharp (5-function in Eq. ( 60 ) , as well as in the analogous 
equation for r+, by a broader probability distribution 



also be broadened broadened by convolution with pi^t in 
both the a and t arguments. 

Finally, the calculations outlined in this paper do more 
than answer a question of principle. As previously men- 
tioned, they serve as a toy model for resolving the chal- 
lenges inherent in the problem of observables in quantum 
gravity. So the lessons learned here, may be applicable to 
a situation like early Universe cosmology, which is a more 
likely source of physically measurable quantum gravita- 
tional effects ^ . 



VII. SKETCH OF CALCULATION IN 
QUANTUM LINEARIZED GRAVITY 

The point of the preceding section was to motivate 
that the output of an explicit quantum calculation should 
be a probability distribution Psit)^ which should be in- 
terpreted as the spectral density of the time register T 
with respect to the quantum gravitational vacuum (pro- 
jected onto the s-eigensubspace of the time register S). 
Phenomenologically, dropping the subscript s since no 
other probability distribution would be considered from 
now on, p(t) should be interpreted as the statistical dis- 
tribution of measurement outcomes for an ensemble of 
repeated measurements of the emission time f (s) (repro- 
ducing the geometry of the experiment for each repeti- 
tion) . Again, motivated by the discussion of the preced- 
ing section, we propose that, within the linearized grav- 
ity approximation and keeping all available parameters 
tuned to minimize all influences on the measurement of 
emission time other than the effects of the gravitational 
field, the role of the observable T should be played by 
the linearized expression ( [SO] ) with the classical graviton 
field everywhere replaced by the quantized graviton field 
(with one caveat to be discussed below) and the role of 
the quantum gravitational vacuum should be played by 
the Poincare-invariant Fock vacuum of the graviton field. 
Within this proposal, the probability distribution p{t) 
can be computed explicitly. We leave the details of this 
calculation to be presented elsewhere [35] and only dis- 
cuss some general aspects of it that can be deduced from 
dimensional analysis and the nature of perturbative cal- 
culations. 



A. Gaussian spectral density 

In linearized gravity on Minkowski space, we interpret 
the Poincare-invariant Fock vacuum |0) as the quantum 
gravitational vacuum and the operator f(s) as the quan- 
tum emission time. From the preceding discussion, our 
goal is to evaluate the spectral density p{t) of |0) with 
respect to t(s). This simplified problem has an explicit 
solution. A linear field theory is essentially a collection 
of harmonic oscillators and, by construction, the Fock 
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vacuum is a Gaussian state (with zero mean) with re- 
spect to the osciUator variables. In other words, the Fock 
vacuum is also Gaussian with respect to any observable 
linear in the graviton field h{x), such as f (s). Therefore, 
the sought probability distribution p{t) is Gaussian. It is 
fully determined by its mean, which is just the classical 
Minkowski space expression (20), and its variance, which 
can be obtained from the expectation value (0|f(s)^|0). 

It is clear that the calculation of the probability dis- 
tribution p(t) is reduced to evaluating a single vacuum 
expectation value given above. Recall that the emission 
time is invariant with respect to gauge transformations 
that fix the synchronization point O and the lab tetrad 
frame at it. However, due to the Poincare invariance 
of the Fock vacuum, the expectation value (0|f(s)^|0), 
which combines the observable and the state, is actually 
invariant under arbitrary gauge transformations and no 
longer depends on the special choice of synchronization 
point or lab tetrad frame. 




T-cl(s) 



FIG. 6. Spectral density of the linearized quantum gravita- 
tional vacuum with respect to the emission time t{s). In a 
linear field theory, the spectral density is expected to be Gaus- 
sian, with mean rci(s) and variance ((Ar)^). The mean is the 
emission time in Minkowski space, Eq. (20 1. This probabil- 



ity distribution clearly penetrates the shaded region, which is 
forbidden by classical causal inequalities. 



B. Causal inequalities and perturbation theory 



In Sec. |IV[ we found that the time delay and emission 
time observables obey some causal inequalities. The va- 
lidity of these inequalities relies mainly on the Lorentzian 
character of the metric tensor and the geodesic charac- 
ter of inertial motion. Therefore, classically, violations of 
these inequalities would be evidence of superluminal sig- 
nal propagation or violation of the equivalence principle. 
Thus we can naturally take the following objective cri- 
terion for the presence of causality violation in quantum 
theory: violation of causal inequalities by the spectral 
density of the quantum gravitational with respect to the 
time delay observable. 

Quantum theory is famous for tunneling phenomena. 
For example, a quantum state may be such that a mea- 
surement may find a particle (though likely with only 
small probability) in a region that is classically forbid- 
den to it. Similarly, we would like to investigate whether 
the causal inequalities are strictly obeyed in the quantum 
theory or are subject to violations via "quantum tunnel- 
ing." A definite answer to this question would go a long 
way toward informing the debate on whether any quan- 
tum theory of gravity necessarily entails causality viola- 
tions [351 133 ■ While it would be very difficult to settle 
this debate, in large part due to the breadth of the sub- 
ject matter, as stated. However, an explicit example of 
a quantum gravitational theory without causality viola- 
tion would force a weakening of the "necessarily entails" 
clause. Equally, a fairly conservative (no extra matter, 
no extra dimensions, no causality violation in the classi- 
cal limit, though taken only in a linear approximation) 
example of a quantum gravitational model with causal- 
ity violation would strengthen the evidence for the "any" 
clause. 

Unfortunately, as should become immediately obvious, 
the perturbative calculations outlined in this section are 



not conclusive enough to establish whether causality vi- 
olation actually takes place or not. In short, since the 
spectral density is expected to be Gaussian (as discussed 
in SeclYll), 



p{t) - exp - 



{t - Tel 

2(Ar) 



(62) 



with some mean Td and variance (Ar)^, all real values 
of t acquire a non-zero probability of being measured. 
Thus, the causal bounds on t are clearly violated, as il- 
lustrated in Fig. [6j However, the responsibility for this 
violation can be ultimately traced back to the perturba- 
tive approximation rather than to the quantum theory. 
Recall that classically, as discussed in Sec. |IV[ the proofs 
of these causal bounds crucially relied on the Lorentzian 
character of the metric, as well as on the detailed be- 
havior of geodesies in Lorentzian spacetimes. Neither of 
these properties survives in perturbation theory. One can 
find classical field configurations of hij{x) for which the 
linearized classical expression for r(s) violates the causal 
inequalities as well. For these field configurations hij{x) 
would have to be of the same order as the background 
metric 77, which is precisely the regime where perturba- 
tion theory is no longer applicable. 

In conclusion, the perturbatively calculated p(t) may 
be presumed to give accurate results around the interval 
[tci — At, Tci -|- At] , but not for larger or smaller val- 
ues. Unfortunately, the information needed to decide 
whether causal inequalities are actually violated requires 
the knowledge of p{t) precisely in the regions where the 
perturbative approximation is no longer expected to be 
valid. 
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C. Finite measurement resolution 

The detailed calculation of the vacuum fluctuation of 
f(s) immediately presents a problem: it is infinite. This 
infinity can be traced back to the singularity of the two- 
point function 



G{x ~y) = {h{x)h{y)) 



1 



{x - yf 



(63) 



in the coincidence limit x — > y. This infinity has 
a straightforward physical interpretation, which at the 
same time suggests a meaningful regularization of the 
divergence. 

Any realistic measurement of the quantum field h{x) 
is carried out by a detector with finite spatial and tem- 
poral resolution. Thus, no measurement is ever sensitive 
directly to the field evaluated at a single spacetime point 
a;, rather measurements are typically sensitive to smeared 
fields 111139] 



(64) 



where g{z) is a smooth test function peaked in the neigh- 
borhood of 0, and ((• • •)) denotes the smearing with re- 
spect to g{z). The smearing function g{z) may be inter- 
preted as the detector sensitivity profile, which clearly 
depends on how the measurement was carried out. The 
vacuum fluctuation of the smeared field at x is then al- 
ways finite 

{h{xf) = / dzidz2 G{x - zi -f- Z2)g{zi)g(z2) (65) 



h{x) — {{h{x — z))) ~ dz h{x — z)g{z), 



dz G{x — z)g{z) 



' {x - zf 



(66) 
(67) 
(68) 



where g{z) = {g * g){z) is the convolution of g{z) with 
itself, by abuse of notation ((• • •)) also denotes smearing 
with respect to g{z), and ^ is the length scale over which 
g{z) has appreciable support, which is the spatiotempo- 
ral resolution of the detector. Physically, this estimate 
means that the root-mean-square noise in a detector, due 
to quantum fluctuations, grows as inversely proportional 
to its resolution ([38] and Sees. 10.9.1-2 of |40]). Such 
fluctuations are vividly illustrated in the context of quan- 
tum optics in Fig. 2.1 of [¥T] . 

Since we are working with an idealized model of phys- 
ical measurement, it is natural that the quantum fields 
entering into the expression for the emission time t{s) 
should be smeared. Unfortunately, the details of pre- 
cisely how the smearing is to be done are quite compli- 
cated. They in general depend on all the aspects of the 
experiment: the resolutions of the proper time clocks, 
the coupling of the lab and probe centers of mass to the 



gravitational field in geodesic motion, sharpness of the 
signals transmitted by the probe, etc. For the purposes 
of this discussion, we do not need such detailed informa- 
tion, as for simplicity we would only be interested in the 
asymptotic limit of perfect detector resolution /i — >■ 0. 
This limit is obviously divergent, so we can settle for 
the leading term in an expansion in inverse powers of /z. 
Therefore, we simply assume that all occurrences of the 
point field hij(x) are replaced by the smeared field h{x)^ 
Eq. ( [64| ). That is, the smearing function g{z) is the same 
everywhere, independent of x. The only thing we assume 
about g{z) is that it is regular enough to render the vac- 
uum fluctuation of f(s) finite and that it is peaked only 
at the origin, with appreciable support over a region of 
size /i, so that we can estimate its moments as 



(69) 



where z*^ represents any homogeneous expression of or- 
der k in the components of z. It is worth noting that, 
as stated, this smearing convention breaks background 
Lorentz invariance. This is clearly unphysical. Neverthe- 
less, we make this assumption in the current and some 
future calculations for the purposes of working out their 
general structure. A more physical smearing convention 
should be re-examined in the future alongside with more 
realistic models of lab, probe, and signal subsystems. 

It is worth noting at this point that the works of Ford 
et al. [12 [TMTT] took a completely different approach to 
the regularization of divergences arising from the singu- 
larities of the graviton two-point function. In particular, 
they treated several scenarios that produced fiuctuations 
different from Minkowski space (finite temperature state, 
squeezed vacuum, extra compactified dimensions), which 
were regularized by subtracting the divergent Minkowski, 
Poincare-invariant vacuum result. Thus these previous 
calculations computed the deviation of the quantum fiuc- 
tuations from that of Minkowski space, but did not di- 
rectly address Minkowski space results themselves, unlike 
we do in this work. 



D. Dimensional analysis 

Looking at the structure of the explicit expression for 
the linearized correction rTc\{s) to the emission t(s), 
Eq. (50), it is fairly obvious that a detailed calculation 



of the variance 



of the smeared correction f, where 



each occurrence of the classical field h{x) is simply re- 
placed by the smeared quantum field h(x), will be quite 
involved. The expression for r, whose structure is illus- 



trated at the end of Sec. ( V B 2 ) , contains on the order of 
10 terms. Therefore, the number of terms in will be of 
order 100. Each of these terms consists of two (possibly 
iterated) integrals over spacetime segments over (possi- 
bly iterated) derivatives of the smeared ((G(z))) graviton 
two-point function. The total number of nested integra- 
tions for each term is five (5), which includes two (2) 
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from the spacetime segment and three (3) from smear- 
ing. Using symmetry, one or two integrations may be 
made trivial. However, it is unavoidable that each of the 
order 100 is a high-dimensional integral. Moreover, the 
integrands are distributions, rather than continuous func- 
tions, whose singularities are ultimately traceable to the 
light-cone and coincidence singularities of the graviton 
two-point functions. The high dimensionality of the in- 
tegrals and the distributional character of the integrands 
makes it very difficult to treat them numerically. On the 
other hand, the integrands of these order 100 terms may 
have many different algebraic structures, preventing the 
evaluation of a single master analytical expression that 
could be uniformly applied to all of them. Splitting each 
term into simpler pieces and considering all possible cases 
of algebraic structures easily leads to thousands of indi- 
vidual integrals to be evaluated analytically. There is 
little choice but to resort to hybrid numerical- analytical 
calculations automated using computer algebra software. 
These detailed calculations are in progress and their re- 
sults will be reported elsewhere [32]. In the rest of this 
section we concentrate on some intermediate, qualitative 
results that may be obtained by straightforward dimen- 
sional analysis. 

Taking dimensionful constants into account, and keep- 
ing in mind that the field h{x) is itself dimensionless, the 
unsmeared graviton two-point functions has the form 



{h{x)h{y)) ^ G{x - y) 



{x - yf 



(70) 



where the denominator of the last expression is the space- 
time interval squared, (x — y)'^ = rjij {x — yYix — y)^ , and 
the numerator is the Planck length squared, — Gh/c" . 
What is important here is that G{z)/£p ^ 1/z^ is a 
homogeneous function of z of degree —2 and hence of 
length dimension [G{z)/£p] = —2. The scales ^ and £p 
and the components of z itself all have length dimension 
[z] = [fj] = [ip] = 1. On the other hand, a derivative with 
respect to z has length dimension [V] = — 1. Generically 
it has the effect Vz" ~ z""^. Using the convention from 
the lAppendix} the spacetime segment integrals are all 
affinely parametrized from to 1 and hence are dimen- 
sionless, [J-^] — 0. On the other hand, integration over a 
spacetime segment has the generic effect z" ~ z"+-'^/s, 
where s is the length scale of the segment X, [s] ~ 1, and 
z on the right hand side corresponds to the coordinates 
of the segment's end points. 

Without smearing, the expectation value (r^) is infi- 
nite. Smearing introduces a regulating length scale fi, the 
detector resolution. Therefore, the smeared expectation 
value (f^) should diverge as /i 0. The details of the 
approach of /i to in general depend on the details of 
the smearing functions. Fortunately, a kind of universal- 
ity among all well-behaved localized smearing functions 
can be obtained by concentrating on the leading terms 
in an expansion of the result in inverse powers of fi. 

From the structure of the explicit expression for r, 
keeping in mind that derivatives worsen singularities 



while integrals improve them, the most singular contri- 
bution should come from the terms with the greatest 
number of derivatives and the least number of integrals. 
Namely, r ^ rx ^^Vh, where rx is some tensorial coefR- 
cient dependent on the geometry of the segment X. Note 
that, since both r and h are dimensionless, the tensorial 
coefficient rx must have length dimension [rx] = 1 and 
be of order s in magnitude, due to the standard affine 
parametrization of the integral over X. In fact, it should 
be of size s, which is the length scale of the spacetime seg- 
ment X. The leading-order contribution to the smeared 
variance of r can then be estimated as follows: 




s'{{ / V-'G(j))) 
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(71) 
(72) 
(73) 
(74) 
(75) 



Detailed calculations show that many terms do have 
this scaling behavior, but also that terms of the form 
(£p//z^) log(/x/s) and (fp//x^)(s//z) show up at intermedi- 
ate stages as well. While the appearance of logarithmic 
scaling is not unusual in quantum calculations, the last 
term is somewhat surprising and, if uncanceled in the fi- 
nal result, may cast serious doubt on the validity of the 
linearized approximation in the regimes of very large s//i 
ratios. This ratio corresponds to that of the spatial and 
temporal extent of the experiment to the resolution of 
the detectors involved. 



From Eq. ( 50 ) , the perturbative correction to the emis- 
sion time r(s) and the time delay (5t(s) scale like sr, and 
so the quantum variances of t{s) and 5f(s) should scale 
like s^(f^), since (f) = 0. From this and the possible 
leading-order contributions to the smeared variance of r 
we can deduce the root-mean-square size of fluctuations 
expected in observations of the time delay due to the fluc- 
tuations of the quantum gravitational vacuum shown in 
Table [T] Let us contrast two possible experimental con- 
texts. In the laboratory context, the spatiotemporal ex- 
tent of the experiment (with time-length conversion via 
the speed of light) is expected to be s ~ 1 m ~ 10~^ s, 
while in the cosmological one s 1 Mpc. Recall that a 
megaparsec is 1 Mpc ^ 10^^ m ^ 10^^ s. For the de- 
tector resolution scale, we select /i ~ 1 nm ~ lO"'^® s. 
This is of the order of the wavelength of X-rays, which 
are consistently available in both contexts. The Planck 



scale as usual is 



10 



-35 



10- 



-44 



All of the above estimates, except one, are well below 
the sensitivity or noise thresholds of the current state 
of the art of experimental and observational technology. 
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TABLE I. Estimates for the root-mean-square size of quan- 
tum fluctuations in observations of the time delay for different 
possible leading-order behaviors in — ^ 0. 



Context 


s 




sip 
M 


^l0£ 






laboratory 


1 m 


1 nm 


10"^^ s 


10" 




10-3" s 


cosmological 


1 Mpc 


1 nm 


10"^^ s 


10" 


" s 


10^ s 



So it is not surprising that kind of effect has yet to be 
observed. Clearly, if the largest of the above estimates 
were correct, we would have observed this effect long ago 
due to the very large fluctuations in the arrival times of 
high frequency photons from distant galaxies. Of course, 
since that result is only preliminary and comes from the 
least understood part of intermediate calculations, it has 
to be taken with a grain of salt. But it does highlight the 
fact that the linearized approximation employed in the 
calculations described above may not be valid on large 
timescales. This is not an unusual feature of perturbation 
theory. For example, it was noticed long ago in celestial 
mechanics that there exist perturbative terms that scale 
with positive powers of time, so-called secular terms, in 
otherwise non-perturbatively stable systems [42]. This 
remark also offers some hope that if, in fact, the pertur- 
bation expansion in our calculations breaks down on large 
time scales that this problem could be repaired using the 
methods already developed for dealing with secular per- 
turbative terms in celestial mechanics or other fields. 



VIII. DISCUSSION 

We have operationally defined a particular physical ob- 
servable, the time delay St{s) [as well as the related emis- 
sion time t(s)], and have provided both exact, implicit 
and approximate, explicit mathematical models for it. 
The time delay satisfies two important inequalities (stem- 
ming from the maximality of light speed and from local 
geodesic extremality) directly related to the causal struc- 
ture of classical Lorentzian spacetimes. Thus, it is sensi- 
tive to the causal structure of classical dynamical gravity. 
Moreover, we have sketched how the same operational 
definition can be used to define a quantum time delay 
observable and how to compute its variance due to quan- 
tum fluctuations of the quantum gravitational vacuum, 
in linearized gravity, given the usual Fock quantization 
of the graviton field. 

This work opens up many potential lines of investiga- 
tion. Foremost among them, is the completion of the 
detailed calculation of the variance of the time delay due 
to the quantum fluctuations of the quantum gravitational 
vacuum. That work is in progress and will be reported 
on elsewhere [32]. 

An important issue that needs to be explored is the 
detailed construction of a quantum model of the mea- 



surement apparatus sketched in Sec. VI This model 



should take into account the quantum dynamics of the 
center of mass motions of the probe and laboratory, a 
more detailed representation of the time stamped sig- 
nal transmitted by the probe, and of the weak measure- 
ments of the relevant clock and signal systems. Some 
existing literature may be helpful in refining these mod- 
els [Ml 1301 131 133] ■ 

The triangular geometry of the time delay experiment 
is one of the simplest possible. However, there is no con- 
ceptual obstacle to generalizing the same methodology to 
more complex geometries, including piecewise geodesic 
motion with more components and even accelerated mo- 
tion. It is also natural to capture other effects of the 
fluctuating gravitational field on the signal, such as an- 
gular blurring and other image effects at the reception 
of the signal by the lab. These effects were previously 
considered in [17], though with caveats similar to those 
given in the Introduction while discussing |12| . 

It is clear that a whole class of physical observables 
of manageable mathematical complexity and with clear 
physical interpretation can be constructed using the same 
methodology. This class can be aptly named astromet- 
ric observables or quantum astrometric observables, when 
referring to them in the quantum context. 

Yet another important generalization is to background 
geometries other than Minkowski space. Cosmological 
and black hole backgrounds are of particular importance. 
For instance, a similar calculation could model the fluctu- 
ation in the arrival time of photons from distant galaxies 
due to the intrinsic quantum fluctuation in the cosmologi- 
cal quantum state of the graviton field. Such fluctuations 
would contribute to the spread of the arrival times of pho- 
tons from distant 7-ray bursts [13]. Undoubtedly, the 
final observational data compounds many effects, includ- 
ing the likely more dominant astrophysical ones and those 
due to in transit scattering. However, a thorough under- 
standing of quantum fluctuations in astrometric observ- 
ables in linearized gravity (or related approximations) is 
necessary before the observational data could be used 
to infer the existence of exotic effects like violation of 
local Lorentz invariance, spacetime discreteness or gran- 
ularity, modified dispersion relations, etc. [I51 - H5] . since 
the model of quantum gravity considered in the present 
calculation exhibits none of these features. Also, the be- 
havior of light signals and inertial or accelerated probes 
in the vicinity of a black hole can be used to give an oper- 
ational meaning to the location of its horizon. The fluc- 
tuations of some quantum astrometric observables could 
then be used to unambiguously study the inferred quan- 
tum fluctuations of the black hole horizon. 

A limitation of the proposed method of calculating 
the quantum vacuum fluctuation of the time delay (or 
any other astrometric observable) in quantum linearized 
gravity is the inability of perturbation theory to ad- 
dress questions involving strong fields, like the question 
of whether the quantum theory respects or violates the 
causal inequalities discussed in Sec. |1V[ Unfortunately, 
in the physically relevant case of four-dimensional space- 
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time, the only effective calculational tool we have is per- 
turbative quantum field theory. Perhaps an improved 
perturbation theory in the spirit of the Magnus expan- 
sion [46' can be used to keep the signature of the metric 
tensor Lorentzian while still using perturbative methods, 
so that the causal inequalities are not immediately vio- 
lated already at the classical level. On the other hand, 
the time delay and astrometric observables in general can 
be defined equally well in any spacetime dimension. This 
opens up the possibility of adapting the quantum calcula- 
tion to the two- and three-dimensional versions of general 
relativity, which can be solved exactly. A family of clas- 
sical observables of 3-dimensional gravity that could be 
said to fall into the astrometric category have been iden- 
tified and expressed in variables that are appropriate for 
treatment in the quantum theory in |47j. The quantum 
calculations have yet to be carried out. 

Finally, since astrometric observables are defined in a 
way independent of the underlying model of quantum 
gravity, their behavior could in principle be studied in 
any of the popular (or even not so popular) proposed 
theories of quantum gravity. It is often the case that it is 
difficult to compare calculations between these different 
theories, due to the very different underlying mathemat- 
ical frameworks. It would be very interesting to see if 
quantum astrometric observables can serve as a bench- 
mark suite to compare the predictions of each of these 
theories on equal footing. 
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Appendix: Perturbative solution of geodesic and 
parallel transport equations 



(Eq. (26)) as el_ = T^x{, the geodesic and parallel trans- 



Let ef be a tetrad field, as described in Sec. VB 1 
Let j(t) be a parametrized spacetime curve and 
a = 0, 1, 2, 3, an orthonormal tetrad along it. Its compo- 
nents (t) in the basis of the spacetime tetrad are given 
by Vait) = ''^ ait) hit))- The pair (7,^°) is a geodesic 
with a parallel-transported orthonormal frame on it if it 
satisfies the following conditions 



at)" = voit)^ 

7(t)'^V,<(t) = 0. 



(A.l) 
(A.2) 



When the spacetime dual tetrad field is expressed in 
terms of a reference inertial coordinate dual tetrad xi 



port equations are expressed in tetrad components as fol- 
lows 



(A.3) 
(A.4) 



where r/kk'^i'^ j = ^ikj = '^i[kj] are the Ricci rotation co- 
efficients (Sec 3.4b of [IS])- The Ricci rotation coefficients 
can be computed in terms of the transformation matrix 
Tj. Below, da — x\di denotes the coordinate derivative, 
F^^ the usual Christoffel tensor, encoding the difference 
between and da, and T^ab = gcc'^lb- 



(^ikj = EkcefVae^j 



Bkcetdae^j + efcce"r^(,e5 
ekcBi OaBj + eje^ej,l cab- 



(A.5) 
(A.6) 
(A.7) 



Each term on the right hand side is evaluated separately 
below and expressed in terms of a single quantity aikj ■ 

a,k, ^ ff {d,,T^,7^ik)ff (A.8) 

ekce^dae'j = ff {doff )T{,7jikX^x^,x'^, (A.9) 

^ff{d,,T^)T{,f^ik (A. 10) 

= -ffffideTf)fi:,Tl,rjik (A.ll) 



-f: id^,T'r]ik)f' ^-a,kj, 



(A.12) 



e-e'^elTcab = ^e^ey^idagbc + db9ac - d^gab) (A.13) 



' rpj rpk' 



X [d.i>gj'k' + dj'gi-k- - dk-gi-y 



(A.14) 



1 



-nf^^n [deiT^nijKT^) (A.15) 
+ dy iTl^iKT^ ) - dk' iTlvijTf,)] 



[aikj + atjk + ajik 

^" ^jki ^kij ^fcjz] 1 
1 



(A.16) 



^jki ^kij ^kji] 

1 , 

+ ajki ~ akij - ttkji] 
= -'^0![ikj] + 2a[j|i|fe] 



(A.18) 

(A.19) 
(A.20) 



The alternative expressions for ujikj in terms of a^kj are 
provided for convenience. When T — exp(/i), and h is 
considered to be small, the linear-order expression for a 
in terms of h is 

a^kJ = d^h%k + 0{h^) = d^hkj + 0{h^). (A.21) 
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7(0) 



7(1) 



7(0) 




7o 



FIG. 7. A geodesic 7 is defined by its initial point 7(0) and 
initial tangent vector 7(0). The initial point itself is specified 
as the final point 7(0) = /3(1) of another curve /3 which starts 
at the origin. The initial tangent vector can then be specified 
by its inverse image 70 G TqM under parallel transport along 
p. 



The geodesic (A.3) and parallel transport (A.4) equa- 



tions can be jointly transformed into a system of integral 
equations 



(A.22) 



t;^(i)-Texp 



exp(p^(t))f<(0), 



viiO), (A.23) 
(A.24) 



where Texp(---) denotes the time-ordered exponential 
and the parallel propagator exp{p^{t))j is defined implic- 
itly by the last equation. For brevity, we also use the 
notation = ^^(l). In this form, the solution can be 
directly expanded to any desired order in 0{h). The so- 
lutions are parametrized by the initial data 7(0)* and 
<(0), withr(0)=«S(0). 

The initial data are specified as described in Sec jVB l| 
Namely, given a curve f3 starting at the origin, /3{0) — O, 
we have /3'(1) = 7'(0) and vl^{0) = exp(p^)*L^w^^^, for 

some vectors Vq ^e^ 
O. 



a _ t3 



^k^o a^j ^^'^ tangent space at 



Suppose that at zeroth order we are given /3{t) = 
/3o(t) + /3i(i) + 0(/i'),7(t) =7o(i) + 0(^),and exp(p)^^ = 
6 J +0{h). To linear order, the parallel propagator is ex- 
panded as 



(^^7o(i))fej-%.(^7o(*))} 

dt'fo{t'M-,o{t')Uj, (A.27) 



(A.25) 
(A.26) 



H,,il). (A.28) 

Recall that L — exp{ho)- The tetrad components of the 
parallel-transported vector at 7(1) are then 



<(1) = exp(p^)} exp(p/3)^ exp(/io); Wo,a 



{he 
0{h^). 



{Hp„ + H. 



7o^i^O,a 



(A.29) 
(A.30) 



For the coordinates of the geodesic curve, the linear-order 
solution is [recall that T = exp(— /i) and L — exp(/io)] 



y (1) = + / dt' eM~Klit')))] exp(P7(i'))i exp(p^)f eMhotn^'Sfl (A.31) 

JO 

= /3^(1) -f- -t- /' dt' [-h){^o{t')) + {H,,{t% + {Hp,)) + {ho)]] v^o + 0{h') (A.32) 
Jo 

= Pl{l) + P\{1) - f dt'h){^/o{t'))vl,^, ~ C dt: f dt"7^-(i")'^(7o(i"))fe^-«'o.o (A.33) 
Jo Jo JO 

+ {H0oyA,o + iho)yo,o + O{h% 

= /3^(1) + /3l{l) + {J^„,pJ + 0{h'), {AM) 

(^7o,/3o)» = mAJlo-PoY (A.35) 

dt'h,,{Mt'))v'o,o- f d<' / dt''^^{t'')uJ{Jo{n)k^Jv'o,0 + iHp„)^,vl,^, + {ho)^,viy,. {AM) 
JO JO 



For the purposes of this paper, all the zeroth-order curves /3o(i) or 70 (t) are piecewise straight line segments in 
Minkowski space. Given a straight line segment X{t), we denote its point of origin and end point by Xi = X{0) 
and X2 = ^(1) respectively. The standard affine parametrization is X{t) = (1 — t)xi + tx2- The segment's standard 
tangent vector is denoted x = X2 — xi. At zeroth order, the /3-parallel-transported image of at O is just Vqq = x^. 
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It is convenient to use the following notation and identities for integration over line segments: 



An) pi rto ptn-i 

/ dtf= dto Ah--- &tnf{X{t^)) (A.37) 
Jx Jo Jo Jo 

dt^-^f{X{t)), (A.38) 

f fiO) 

dtf= dtf, (A.39) 
Jx Jx 



[f]ll = [ dtx%f = f{x2) - f{xi), (A.40) 
Jx 

/ dtx'dif = .f{xi)+ dtf. (A.41) 

Jx IT'- Jx 



For definiteness, suppose that 70 is a single segment X, whose point of origin is the end point of /3o, a piecewise linear 
path Y, whose segments are indexed by A'' and denoted Yn- Then we can write the expressions for Hx, Hy and Jx,y 
more concretely as 

A'B^{Hx)ij=A'B^ [ dtx''[dkhiij]-djh^ki)+dih^kj)] (A.42) 
Jx 

= A'B^ [h^ijf.l + 2A^'B^^x'' [ dt dih^kj) , (A.43) 

Jx 

{HY)ii =Y.{Hy,)ii, (A.44) 

N 

where the expression for ily^ is the same, only with X replaced by Yjv, and 

A\Jx,Y)^ = -A' / dtKjX^ + A' / dtx'^ldkhy.j} - ^,h^k^) + d,h^kj))]x' + A\HY)ijX^ + A\ho)ijX^ (A.45) 
Jx Jx 

= A'x^ {ho)ij - A'x^h^i^} + A'x^ [h^ij^]^' (A.46) 

N JYn 
r-(l) 



- A'x^ I dt + 2A^^x^^x'' ^ dih^kj) + Yl 2^''^^^^ / dtdih^ 
A'x^ [{hoh - hj] (O)] - A"x^ I + 2A^'x^x'' t ^ dt dih^kj) + 2A^'x'^y% [ dt dih^kj) ■ 

Jx Jx »r JYm 



(A.47) 

Given the above explicit formulas for Hx and Jx,y- We can consider how they transform under gauge transfor- 
mations. In the tetrad formalism, gauge transformations are generated by arbitrary local Lorentz transformations 
(tetrad rotations that do not change the metric) and arbitrary spacetime diffeomorphisms. In the linearized tetrad 
formalism, the most general gauge transformation takes the form hij hij +djCi+Dij and {ho)ij ^ — SyCjj +Dij, 
where Dji = —Dij. Under such a variation we have the following identities: 

Shij=djCi + Dij, (A.48) 
6{ho)ij=2dyC^+Dij, (A.49) 

A'Bi {SHx)ij = A'Bi [^(7^] + A,]^? + ^''-B^l / dt [x^dkdkCj + x^didjCk] (A.50) 

Jx 

= A'Bi[d^C,^ + D,,]ll + A^^B'\d,C,]ll (A-51) 

= A'B^[D,^]%1, (A.52) 

A'{5Jx,Y)i = A'x^[2d[jCi]{0) - - B„(0)] - ^A^^ [ dt [d,C\ + d,Cj] (A.53) 

^ Jx 

+ A^'x^^xf" / dt [didkCj + didjCk] + V A^'x^Vn / dt [didkCj + didjCk] 

Jx jy JYm 
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X 



A^x^dyQ^iO) - B,,iO)] - Ai'x^^d,C,{0) - A\C,f"^ 
A^x\2dyfi,^{0)-B,,{p)\-A\C,\%\. 



(A.54) 



(A.55) 
(A.56) 



Suppose that we have a sequence of segments (X) that 
starts at O and forms a closed loop, = 0. We can 

naturally form pairs [X, F) where the segments Ijv of Y 
consist of elements of [X) that precede X. Then we can 
define 



X 



iA.57) 



{J(X))i = ^{Jx,Y)t- 



X 



(A.58) 



From the above transformation properties, it is clear that 
-ff(x) and J(x) are invariant under all gauge transfor- 
mations, A^B^SH(^x))ij = and A^{SJ(^x)h = 0. The 
terms that depend only on the values of C and D at O 
cancel because their sum is proportional to ^'^^ 
the remaining terms cancel because they form a cyclic 
telescoping sum. These invariant quantities are actually 
used in Sec. |VB 2 in the expression for the correction to 
the time delay in approximately Minkowski spacetime. 
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